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Abstract 

We introduce a one-parameter generalized oscillator algebra A K (that covers the case 
of the harmonic oscillator algebra) and discuss its finite- and infinite-dimensional repre- 
sentations according to the sign of the parameter k. We define an (Hamiltonian) operator 
associated with A K and examine the degeneracies of its spectrum. For the finite (when 
k < 0) and the infinite (when k > 0) representations of A K , we construct the associated 
phase operators and build temporally stable phase states as eigenstates of the phase opera- 
tors. To overcome the difficulties related to the phase operator in the infinite-dimensional 
case and to avoid the degeneracy problem for the finite-dimensional case, we introduce 
a truncation procedure which generalizes the one used by Pegg and Barnett for the har- 
monic oscillator. This yields a truncated generalized oscillator algebra A K)S , where s 
denotes the truncation order. We construct two types of temporally stable states for A KiS 
(as eigenstates of a phase operator and as eigenstates of a polynomial in the generators 
of A K , S ). Two applications are considered in this article. The first concerns physical re- 
alizations of A K and A Ka in the context of one-dimensional quantum systems with finite 
(Morse system) or infinite (Poschl-Teller system) discrete spectra. The second deals with 
mutually unbiased bases used in quantum information. 



1 Introduction 



It is well known that the usual model for the quantized single modes of the electromag- 
netic field is the harmonic oscillator with an infinity of states. The infinite-dimensional 
character of the representation space of the corresponding oscillator algebra constitutes 
a drawback to define a phase operator in a consistent way Q]|-ll3]|. In order to get rid 
of this difficulty, Pegg and Barnett suggested to truncate to some finite (but arbitrarily 
large) order the infinite-dimensional representation space of the oscillator algebra @J. 
Their approach also provided a valid way for calculating the so-called phase states (the 
eigenvectors of the phase operator). In the same vein, Vourdas proposed a definition of 
a phase operator for su(2) and calculated its eigenstates without a truncation procedure 
since su(2) admits finite-dimensional unitary irreducible representations [5]. He also con- 
structed a phase operator and its eigenstates for su(l, 1), without a truncation procedure 
although su(l, 1) admits infinite-dimensional unitary irreducible representations flU. 

The main aim of the present work is to develop a method to build unitary phase op- 
erate Jj] and temporally stable phase states for some exactly solvable quantum systems. 
Various algebraic structures were used to construct (temporaly stable or not) coherent 
states in connection with some quantum systems [|6l- [fT0l . The construction of temporally 
stable phase states to be developed in this work is based on a generalized oscillator alge- 
bra which takes its root in [fTTl[T2l . This algebra was introduced to construct isospectral 
shape invariant potentials in the framework of fractional supersymmetry. 

A second facet of this work is to show that the obtained temporally stable phase states 
can be used to generate mutually unbiased bases (MUBs). Such bases are of considerable 
interest in quantum information and were recently investigated from an angular momen- 
tum approach lfT3l[T4l . It is not the purpose of this paper to deal with unsolved problems 
concerning MUBs but to give a way to construct MUBs from temporally stable states 
associated with some exactly solvable systems. 

The paper is organized as follows. Section 2 is devoted to the generalized oscillator 
algebra A K . Temporally stable phase states associated with A K are studied in section 
3. Section 4 deals with the truncated oscillator algebra A K , S and the correponding phase 
states. As a first application, the derivation of MUBs from phase states is developed in 
section 5. A second application is made in section 6 to some exactly solvable quantum 
systems. 

'We deal here with unitary rather than Hermitian phase operators. The two kinds of operators are related 
via an exponentiation trick. 



2 



The notations are standard. Let us simply mention that: 5 a ,b stands for the Kronecker 
symbol of a and b, I for the identity operator, A* for the adjoint of the operator A, and 
[A, B] and {A, B} for respectively the commutator and the anticommutator of the opera- 
tors A and B. We use a notation of type for a vector in an Hilbert space and we denote 
(4>\ip) and respectively the inner and outer products of the vectors and \<p). 

2 Generalized oscillator algebra 

2.1 The algebra A K 

Let A K be the algebra spanned by the three linear operators a~, a + and N satisfying the 
following relations 

[a",a+] = I + 2kN [N,a ± }=±a ± (a - )* = a + iV f = N, (1) 

where k is a real parameter. Note that, for k — 0, the algebra Aq is nothing but the 
usual harmonic oscillator algebra. The operators a~, a + and iV in (Q]) generalize the 
annihilation, creation and number operators used for the harmonic oscillator. Therefore, 
the algebra A K shall be called generalized oscillator algebra. This algebra turns out to be 
a particular case of the generalized Weyl-Heisenberg algebra Wk introduced in |fTTl [T2l| 
and not to be confused with the Lie algebra of the Heisenberg-Weyl group HW(M) used 
in quantum information [[141 . In fact, A K is identical to Wk with 

k = l f (N) = aN + b -±=X ± = a ± « = ^p (2) 

where the operators fo(N) and X±, and the parameters k, a and b are defined in |fT2|. It 
should be noted that the C*A-extended oscillator algebra worked out in lfT51 is a particular 
case of Wk (for A = k). 

2.2 Hilbertian representation of A K 

We denote by T K the finite- or infinite-dimensional Hilbert space on which the operators 
ar , a + and are defined. Let 

{\n) :n=0,l,...,d(«)} (3) 

(with din) finite or infinite) be an orthonormal basis, with respect to the inner product 
(n\n'} = 8 n ^ n i, of the space T K . It is easy to check that the actions 

a + \n) = v^VTT)e- i[F(n+1) - F(n) ^|n+ 1), 
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a -\n) = ^Finje 1 ^-^- 1 ^ - 1), (4) 
a~|0) = N\n)=n\n) 

provide an Hilbertian representation of the algebra A K defined by (CO)- In equation ©, the 
real parameter ip is arbitrary and the positively-valued function F : N — > R + satisfies the 
recurrence relation 

F(n + 1) - F(n) = l + 2nn F(0) = 0. (5) 
The iteration of © yields 

F(n) =n[l + K(n- 1)], (6) 

which is linear in n only for k = 0. Since F(n) G R+, we must have the following 
condition 

1 + K (n - 1) > (7) 

for n > 0. The condition © determines the value of d(n) and then the dimension of F K . 
The finiteness or infiniteness of T K depends on the sign of the parameter k. For k > 0, the 
space JF K is infinite-dimensional. In fact, for k = 0, the space JFq coincides with the usual 
Hilbert-Foch space for the harmonic oscillator. For k < 0, there exists a finite number of 
states satisfying the condition ©. As a matter of fact, for k < 0, n can take the values 

n = 0,l,...,S(--) = d-l, (8) 

where stands for the integer part of x. The finiteness of the space JF K induces 

properties of the operators a~ and a + which differ from those corresponding to an infinite- 
dimensional space. In particular, the trace of any commutator in the finite-dimensional 
space must be zero. This implies that the parameter k is related to the dimension d of the 
space T K by 

d = 1 - -. (9) 

K 

Equation © requires that — 1/k be a positive integer. In the following, we shall assume 
that e W when k < 0. 



4 



2.3 A generalized oscillator Hamiltonian 

We are now in a position to define an operator which generalizes (up to an additive con- 
stant) the Hamiltonian a + a~ + 1/2 for the one-dimensional harmonic oscillator. Starting 
from 

a + a~\n) = F(n)\n) F(N) = a + a~, (10) 

we refer F(N) to as an Hamiltonian associated with the generalized oscillator algebra 
A K . The eigenvalue equation 

F(N)\n) = n[l + k(u - l)]\n) (11) 

gives the energies © of a quantum dynamical system described by the Hamiltonian op- 
erator F(N). Let us discuss the degeneracies of the levels F{n) given by ©. 

(i) In the case k > 0, the spectrum of F(N) is nondegenerate. 

(ii) In the case k < 0, the eigenvalues of F(N) can be rewritten as 

d — n 

F(n) = nj-j, (12) 

so that 

F(n) = F(d-n) n = 1, 2, . . . , d - 1. (13) 

Thus, for d even the levels are doublets except the fundamental level n = and the level 
n = d/2 which are nondegenerate. For d odd the levels are two-fold degenerate except 
the fundamental level n = which is a singlet. 

In both cases (k > and k < 0), we note that the Perron-Frobenius theorem lfT6l is 
satisfied, namely, the fundamental level is nondegenerate. 

It is known that one-dimensional quantum dynamical systems (on the real line) corre- 
spond to nondegenerate spectra. Therefore, the representation obtained for A K with k < 
cannot be used to describe a particle evolving in some nonrelativistic potential on the real 
line. However, a modification of the generalized oscillator algebra A K can be achieved 
in orded to avoid the degeneracies of F(N). This will be done in section 4 by means 
of a truncation procedure which will prove also useful in the case k > to define in a 
consistent way the phase operator for some exactly solvable systems. 

3 Temporally stable phase states for A K 

We shall treat separately the cases k > and k < associated with the infinite- and the 
finite-dimensional representation of the generalized oscillator algebra A K , respectively. 
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3.1 The infinite-dimensional case 

In the case k > 0, we decompose a~ and a + as 

a~ = E^^/fiN) a + = y/F{N) (E 00 ) t , (14) 

where 

oo 

Eqo := J2e i[F{n+1) - F{n)]<p \n)(n + l\. (15) 

n=0 

It is important to emphasize that 

oo oo 

£oo (^oo) 1 " = X] l n >H = 1 (E^E^ = \n){n\ =1- |0)(0|, (16) 

n=0 n=l 

a result which means that is not a unitary operator. 

To find the phase states corresponding to k > 0, let us consider the eigenvalue equa- 
tion 

E OB \z) = z\z) zeC. (17) 
By expanding the vector \z) of T K as 

oo 

\z) = J2C n z n \n), (18) 

n=0 

it is easy to see that the complex coefficients C n satisfy the relation 

C n+l = e -i\F{n+l)-F{n)W Cn fl E N. (19) 

It follows that 

C n = e~ iFin)lp Co neN, (20) 

where the coefficient C can be determined from the normalization condition of the states 
\z). As a result, we can take (up to a phase factor) 

oo 



I*) = yjl - \z\ 2 ^z n e~ lF{n)lp \n) (21) 
on the domain {z e C, \z\ < 1}. 



n=0 
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Following the method developed in [fTVll for the Lie algebra su(l, 1), we define the 
states \9, <p) by 

\0,<p):= lim JL^ k), (22) 



where 6* G [— 7T, +7r] (see also [fT8l where a limit of type z — > e ld =>• |z| — > 1 is used in a 
similar way). We thus get the states 

oo 

\9,(p) = Y,e ine e~ iF(n)v \n). (23) 

n=0 

These states, defined on the unit circle S 1 , turn out to be phase states. Indeed, we have 

E^O,?) =e*\9,<p). (24) 

Hence, the operator is a (nonunitary) phase operator. 
The main properties of the states \9, cp) are the following. 

(i) They are temporally stable in the sense that the relation 

e- iF{N)t \9,ip) = \9,cp + t) (25) 

is satisfied for any value of the real parameter t. This property is due to the presence of 
the parameter (p in the phase operator E^. 

(ii) They are not normalized and not orthogonal. However, for fixed cp, they satisfy 
the closure relation 



1 

2^ 



d9\9,<p)(9,(p\ = I. (26) 



Finally, observe that for <p = the states \9, 0) have the same form than those derived 
in 03 for su{l, 1). 



3.2 The finite-dimensional case 

For k < with — 1/ k E N*, the Hilbert space T K is d-dimensional with d — 1 — 1/ k. The 
action of a~ and a + on T K is given by © supplemented by 

a + |d-l)=0, (27) 

which easily follows from the calculation of (d— l\a~a + \d — 1). 
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Let us look for a decomposition of the creation a + and annihilation a operators 
similar to (fT4)) for the case k > 0. Thus, let us put 



a- = E d ^F\N) ^a + = y/F{N) {E d ) ] . (28) 
The operator E d can be seen to satisfy 

E d \n) = jim-nn-V]<P\ n _ i) ( 29) 

for n — 1, 2, . . . , d — 1. For n = 0, we shall assume that 

£ d |0) = e *[^(o)-^-i)]^| rf _ x) ( 30) 

so that (|29|) is valid modulo d. (Note that, in view of (|28l) . a" |0) = does not imply that 
£J d |0) =0.) It follows that we have 

(E d ) ] \n) = e^K 1 )-^^ + 1), (31) 

where n + 1 should be understood modulo d. As an important result (to be contrasted 
with the situtation where k > 0), the operator is unitary. Therefore, equation (|28l) 
constitutes a polar decomposition of a~ and a + . 

We are now ready to derive the eigenstates of the operator E d . Let us consider the 
eigenvalue equation 

d-l 

S„|z)=z|z) |z) = ^C n z n |n) (32) 

n=0 

with 2 G C. Here again (as in the case k > 0), we obtain a recurrence relation for the 
coefficients C n , viz., 

C n = e-m^-Hn-i^c^ n = l,2,...,d-l (33) 
with the cyclic condition 

C = z d e^ F ^ F{d - l ^C d ^. (34) 

Therefore, we get 

C n = e- iF{n)v C n = 0,1,..., d-l, (35) 
with the discretization condition 

z d = 1. (36) 
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As a consequence, the complex variable z is a root of unity given by 



z = q m m = 0,1,..., d-1, (37) 

where 

q := e 2ni/d (38) 

is reminiscent of the parameter used in the theory of quantum groups. The constant C 
can be calculated from the normalization condition (z\z) = 1 to be 

c ° = 73 <39) 

up to a phase factor. Finally, we arrive at the following eigenstates \z) = \m, ip) of E d 

d-1 

\m,<p) = —y^ e - lF ^q mn \n). (40) 

The states |m, ip), labeled by the parameters m E Z/dZ and ip e R, satisfy 

^ d |m,^) = e l9m |m,v?) ^ m = m— , (41) 

a 

which shows that Ed is indeed a phase operator. In the particular case ip = 0, the states 
|m, 0) are similar to those derived in ifTTl for the Lie algebra su(2). In this case, the 
states \m, 0) correspond to an ordinary discrete Fourier transform of the basis {|n) : n — 
0,1, ... ,d — 1} of the (i-dimensional space T K . 

The phase states \m, <p) have remarkable properties (to be compared to those for the 
states \6, <p) of the case k > 0). 

(i) They are temporally stable under "time evolution". In other words, they satisfy 



e 



- iF{N)t \m,<p) = \m,<p + t). (42) 



for any value of the real parameter t. We note here the major role of the parameter ip in 
ensuing the temporal stability of the states \m,cp). 

(ii) For fixed (p, they satisfy the equiprobability relation 

|(n|m, (p)\ = — -= n, m e Z/dZ. (43) 
yd 

(iii) For fixed if, they satisfy the orthonormality relation 

(m, (p\m', ip) — 6 mm > m,m'£Z/dZ (44) 
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and the closure property 

d-l 

^2\m,<p)(m,<p\= I. (45) 

m=0 

(iv) The overlap between two phase states \mf, if') and \m, if) reads 

1 d-l 

(m, v |m', = 3 E g p(m - mW ' n) , (46) 

n=0 

where 

p(m — m', — n) = — (m — m!)n H (y> — ip')F(n) (47) 

27T 

and g is defined in (1381) . Therefore, the temporally stable phase states are not all orthogo- 
nal. 

4 Truncated generalized oscillator algebra and phase states 

As discussed in section 2, in the case k > the Hilbert space T R associated with A K 
is infinite-dimensional. It is then impossible to define a unitary phase operator (see sec- 
tion 3). On the other hand, in the case k < with — 1/k G N* the space T K is finite- 
dimensional and there is no problem to define a unitary phase operator. However, the spec- 
trum of the Hamiltonian F(N) associated with A K for —1/k E N* exhibits degeneracies. 
Therefore, it is appropriate to truncate the space T K for both k > and k < in order to 
get a subspace F KyS of dimension s with the basis {\n) : n = 0, 1, . . . , s — 1}. For k > 0, 
the truncation is done at s sufficiently large (note that the difference F(n + 1) — F{n) 
between two consecutive states increases with n for k > so that we can ignore, in a 
perturbative scheme, the states with n large). For k < 0, the truncation can be done at 
s — (d + 2) /2 for d even and at s — (d + l)/2 for d odd (with d given by ©) in order to 
avoid the degeneracies of F(N). 

4.1 The truncated algebra A KjS 

Inspired by the work of Pegg and Barnett flU, we define the truncated generalized oscilla- 
tor algebra A K , S through the three linear operators b~, b + and iV satisfying the following 
relations 

[b~,b + ] = I + 2kN-F{s)\s-1)(s-1\ [N,b ± } = ±b ± (b-y = b + iV f = iV.(48) 
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The algebra A K>S generalizes the one introduced by Pegg and Barnett for the harmonic 
oscillator in their discussion of the phase operator for the single modes of the electromag- 
netic field H. Indeed, the algebra Aq >s , for k = 0, is identical to the truncated oscillator 
algebra considered in [4] . 

Following the same approach as in subsection 2.2, we define a s-dimensional repre- 
sentation of A KiS (whatever the sign of k is) via the actions 

b+\n) = y/F(n + l)e- i[F{n+1) - F ( n)]v> \n + 1), 

b-\n) = V 'F{n)e llF{n) - F(n ~ 1)]ip \n - 1), (49) 

&-|0)=0 6 + |s-l) = N\n)=n\n) 

for n = 0, 1, . . . , s — 1. Note that a further condition is necessary here, namely, the upper 
limit condition b + \s — 1) = 0. It can be checked that the recurrence relation © is equally 
valid for A KtS - Therefore, equations © and (PT21) can be applied with n = 0, 1, . . . , s — 1. 

It is interesting to note that the creation and annihilation operators b~ and b + satisfy 
(in the representation under consideration) the nilpotency relations 

(b-y = (b + y = o, (50) 

which are similar (for s = k 6 N \ 0, 1) to those describing the so-called fc-fermions that 
are objects interpolating between fermions (for k = 2) and bosons (for k — > oo) lfT9l . 

4.2 Phase states for A KjS 

For the truncated algebra A K>S (corresponding to d(n) finite or infinite), the analog of the 
phase operator E d is the unitary operator 

s-l 

E s := e iim-F(s-i)]<p\ s _ 1)(Q | + ^[FW-i^-i)]^ _ ^ n |_ (51) 

n=l 

By using the same reasoning as in subsection 3.2, we obtain 

27T 

E s \m, (f) = e %Bm \m,Lp) 9 rn = m — , (52) 

s 

where 

1 

\m^) = J2 e ~ lF{n)iP (ls) mn \n), (53) 
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with m G Z/sZ, <p G R and g s given by 

:= e 2m/s . (54) 

We are thus left with phase states \m, cp) associated with the phase operator E s . These 
states satisfy the same properties as those for E d (see section 3.2) except that d is replaced 
by s in some places. 

4.3 A new type of discrete phase states 

It is well known that, for quantum systems with a finite spectrum (like the Morse system) 
or for Lie algebras with finite-dimensional unitary representations (as for instance su(2)), 
the construction of coherent states cannot be achieved by looking for the eigenstates of an 
annihilation operator jH or of a compact shift operator J6J. 

For the algebra A KtS the difficulty inherent to the finiteness of the representation can 
be overcome as follows. We define the operator 

v -^"' + w=Ty <55) 

where the function E is defined via 

E(0):=1 E(n) := F(1)F(2) . . . F(n) n = 1, 2, . . . , s - 1. (56) 
The operator V s is an idempotent operator of order s since 

{V s ) a = I. (57) 
Let us consider the eigenvalue equation 

s-1 

V s \z) = z\z) \z) = J2°nZ n \n) zeC. (58) 

n=0 

By using (1571) . we obtain that z is discretized as 

z = (q s y fj, g Z/sZ, (59) 

with q s defined by (|54l) . Then, it is a simple matter to calculate the coefficients C n and to 
normalize the \i- and (/^-dependent states \z) = cp). This leads to 



s-1 

\p,<p) = C ^^_(gs)^e-^|n), (60) 

n=0 VH n ) 
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where the normalization factor C is such that (up to a phase factor) 



hi in) 



The states if) are temporally stable and are similar to the coherent states introduced 
by Gazeau and Klauder [ 8 ] except that their labeling includes an integer and they corre- 
spond to the eigenvectors of a polynomial in terms of generalized creation and annhilation 
operators. They satisfy 

(//, <ptf, <p') = Cl £ -^-( qs r^)e- iF ^'-^ (62) 

and 

1 s ~ * ^ 1 

-^\^v){^^\^Cl^—^\n){n\. (63) 

m=0 n=0 ^ ' 

Consequently, they are not orthogonal. 

We close this subsection with a remark concerning the unitary operator 

U s := (q s ) N (64) 

that is a companion of V s in the following sense. This operator satisfies the cyclicity 
relation 

(U s ) s = I. (65) 
Furthermore, we have the s-commutation relation 

V S U S - q s U s V s = 0. (66) 

Equations (1571 ), (|65l) and (|66l ) are necessary conditions for the pair (U s , V s ) be a pair of 
Weyl (see IfPfl ). However, this is not the case because V s is not unitary. 



5 Application to mutually unbiased bases 

As an a priori unexpected connection, the approach in subsection 3.2 and 4.2 for the 
finite-dimensional cases (for A K and A KtS ) can be further developed for deriving MUBs. 
Let us recall that two orthonormal bases {\aa) : a = 0, 1, . . . , d — 1} and {\b/3} : /3 = 
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0, 1, . . . , d — 1} in a ci-dimensional Hilbert space (with an inner product ( | )) are said to 
be mutually unbiased iff 

\(aa\bf3) \ = 5 afi 5 afi + -^=(1 - 5 0>6 ). (67) 

For fixed d, it is known that the number J\f of MUBs is such that N < d + 1 and that the 
limit J\f = d + lis reached when d is the power of a prime number 11201 12T1 . 

5.1 MUBs from phase states for A K 

In order to generate MUBs along the line of the developments of subsection 3.2, let us 
further examine some properties of the phase operator E d for A K with n < 0. This 
operator can be written in a compact form as 

E d = STe^-^-^n - l)(n| (68) 

71=0 

(in this section, the summations on n are understood modulo d). It is easy to check that 

(E d ) d = I, (69) 

so that Ed is idempotent. The operator E d can be decomposed as 

E d = U V V, (70) 

where the operators U v and V are defined by 

d-l 

U v := e mN+i)-F(N) ]lf y \ n - l)(n|. (71) 

ra=0 

The operators and V are unitary and satisfy the pseudo-commutation relation 

U V V = e 2iv/{d - 1} VU v . (72) 
In addition, the operator V satisfies the idempotency relation 

V d = I (73) 
and, when the parameter cp is quantized as 

<p = -n^—-^p P e Z/dZ, (74) 
d 
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we have 

(U v ) d = e i7r{d - 1)p I. (75) 
In view of ( 1741 ), equation (1721 ) can be rewritten as 

VU V = q p U v V (76) 



(see (1381) for the definition of q). For the discrete values of p afforded by (1741) , equation 
(|40l) yields the phase states \m, <p) = \m,p) given by 



l m ' p ) = <f (d ~ n)p/2+m >> p,m G Z/dZ, (77) 



n=0 



which coincides with the vector \aa), with a = p and a = m, obtained in [fT3l in an 
SU(2) approach to MUBs. Alternatively, by putting 

k:=d-n-l \n) = \d-k-l) = \k), (78) 

equation (1771 ) becomes 

d-l 

\m )P ) = J_^2 q (k+W-k-i)p/2-(k+i)m\ k ^ p,m £ Z/dZ, (79) 

which coincides with the vector \aa), with a = p and a = m, derived in |fl4| in an 
angular momentum approach to MUBs. It is to be observed that (|77l) and (|79l) correspond 
to quadratic discrete Fourier transforms. 

To make a further contact with [fT3l[14l . let us note that when ip is discretized accord- 
ing to (1741 ), the inner product (m, ip\m' , ip') = (m,p\m',p') (see equation (|46l ) can be 
rewritten as 

(m,p\m' ,p') — -S(u,v,w) (80) 

CI 

with 

u:=p — p' v := — (p — p')d + 2(W — m) w := d. (81) 

In equation (f8~0l) , the factor S(u, v, w) denotes a generalized quadratic Gauss sum defined 
by [1221 

|io|— 1 

S(u,v,w) := ^ e i7r(ufc2+,,fc)/u ', (82) 

fc=0 

15 



where u, v and w are integers (the nonvanishing of S(u, v, w) requires uw + v even). In 
the special case where d is a prime integer and p' ^ p, the calculation of S(u,v,w) in 
(f8~0l) through the methods developed in ||22l|23) (see also OS) leads to 

\(m,p\m',p')\ = -y=. (83) 

This result shows that the d bases 

B p := {\m,p) : m = 0, l,...,d- 1} p = 0, 1, . . . , d - 1 (84) 

of the d-dimensional space T R , with d given by ©, are mutualy unbiased. On the other 



hand, in view of (1431) . it is clear that any basis B p and the basis 

B d :={\n) :n = 0,l,...,d-l}, (85) 

known as the computational basis in quantum information and quantum computation, are 
mutually unbiased. As a conclusion, for d prime, the d bases B p with p — 0, 1, . . . ,d — 1 
and the computational basis £>d constitute a complete set of d + 1 MUBs. This result, in 
agreement with the one derived in lfT3l[T4l . is the starting point for constructing MUBs in 
power prime dimension. 

5.2 MUBs from phase states for A K , S 

By applying a discretization procedure similar to the one introduced in subsection 5. 1, we 
can construct MUBs from the phase states (1531) for the truncated algebra A K>S with k ^ 0. 
Let us quantize the parameter <p by putting 

2 

m = TT—p p E Z/sZ. (86) 
SK 

Then, equation (|53l leads to the states |m, <£>) = |m,p) given by 



n=0 

where 



1 s_1 

p) = -= V(g s ) n ( 5 - n ) p+nm |n) p, m G Z/sZ, 



(87) 



5 := 1 - -. (88) 

In this subsection, we shall assume that 1/k E Z so that 5 E Z (note that 5 = <i for k < 0). 
The overlap (m, p|m', p') can be written in terms of the generalized quadratic Gauss sum 
(f8~2l) . Indeed, we have 

1 5-1 1 

(m,p|m',p') = - "y( q \n(S-n)( P '- P )+n( m '-m) = _ S f uvw \ (89) 
n=0 
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where 



u := 



2(p-p') 



v := 



25(p — p) + 2(m' — m) 



w := s. 



(90) 



We can proceed as in subsection 5.1 in order to show that the various states \m,p) gener- 
ate, together with the s-dimensional basis {\n) : n = 0, 1, . . . , s — 1}, s + 1 MUBs when 
s is a prime integer. 

6 Application to exactly solvable potentials 

The main goal of this section is to show how the generalized oscillator algebra A K is 
relevant for the study of one-dimensional exactly solvable potentials in the context of 
supersymmetric quantum mechanics and how MUBs can be derived from the temporally 
stable phase states for some quantum mechanical systems. 

6.1 Creation, annihilation and transfer operators 

Ordinary supersymmetric quantum mechanics can be presented in differnt ways (e.g., see 
ll24l - [[30l '). We adopt here the approach according to which a supersymmetric dynamical 
system is defined by a triplet (H,Q + ,Q-) 2 of linear operators acting on a Z 2 -graded 
Hilbert space H and satisfying the following relations 



H = H^ Q- = Q{ Ql = {Q.,Q + } = H [H,Q ± }=0. (91) 



(In this approach, ordinary supersymmetric quantum mechanics is a particular case, cor- 
responding to k = 2, of fractional supersymmetric quantum mechanics of order k dealing 
with triplets (H,Q + ,Q_)k which satisfy relations generalizing (l9TT) and which corre- 
spond to a Zfc grading [12].) The operators Q + and Q_ are the supercharges of the 
one-dimensional system. We suppose that the spectrum of the self-adjoint operator H, 
the supersymmetric Hamiltonian of the system, is discrete. The Hamiltonian H can be 
written 



where H and Hi act on the states |\l/ n , 0) and \^ n , 1) of even and odd grading, respec- 
tively. In other words, the Hilbert space H is decomposed as 



H = H ®Hi Ho := {|*n,0) : n ranging} Hi := {|tf n , 1) : n ranging}, (93) 



H = H + Hi, 



(92) 
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which reflects the Z 2 grading. We shall assume that there is no supersymmetry breaking. 
In this case, the Hamiltonians H and Hi are isospectral except that the ground state of 
H has no supersymmetric partner in the spectrum of Hi. 

By combining the above-mentioned considerations on supersymmetry with the Infeld 
and Hull factorization method [|28l , we can construct creation, annihilation and transfer 
operators for an exactly solvable Hamiltonian in one dimension Il24ll - [|35l . For this pur- 
pose, let us consider a one-dimensional quantum system embedded in a real potential 
v o : x y-¥ v (x). The corresponding Hamiltonian is 

(94) 

Let us suppose that the Hamiltonian H is exactly solvable and admits the discrete spec- 
trum 

e = < ei < e 2 < . . . < e n < e n+ i < . . . , (95) 

with a finite or infinite number of levels. We know that the Hamiltonian H of this system 
can be factorized as |j25l |26l |28ll29l [30l 



+ - + I f d \ _ 1 / d \ 

Hn = x x x := —7= — - — h w ) x := —7= - — h w . (96) 

V2\dx J y/2\dx J 

The superpotential w : x i-> w(x) satisfies the Ricatti equation 

1 / o dw\ 
= n U - -r- • (9V) 



2 V dx J 

Since the ground state energy is assumed to be zero, it is easy to see that the potential v 
and the superpotential w can be expressed in terms of the ground state wavefunction. 

It is important to stress that the operators x + and x~ are not in general creation and 
annihilation operators for Ho ll25l l26l |29l l30l l35l . They are indeed transfer operators 
from the spectrum of H to the one of Hi and vice-versa. To identify them, we start 
by representing the supercharge operators and the supersymmetric Hamiltonian by 2 x 2 
matrices (23 M M ED 



x J ^ + \x + OJ \ H 0/ 

where the operator 

Hi = x -x+ = ~-f-+vi (99) 
2 dx z 
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is the supersymmetric partner of H and corresponds to a new potential V\ : x \-> Vi(x). 
The potential 

is the supersymmetric partner of the potential v . The Hamiltonian Hi is also exactly 
solvable and isospectral to H (except for the ground state). Indeed, 

Ho\^ n ,0} = e°|* n ,0)=J-ff 1 (a;-|* fl ,0)) = e°(a;-|* n ,0)) J (101) 

where e° := e n . Similarly, 

Hx\* n , 1) = e*|tt n , 1) floM*™ 1)) = 4(x + |^ n , 1)). (102) 

(For more details see 11251 l26l l29l and the recent topical review 11301 .1 From equations 
(| 101b and (11021) . it is clear that we can take 

x-\^ n+1 ,0) = v ^Te i ( e " +1 - e "^|^ n ,l) (103) 



a; + |^ n ,l) = v ^T7e- i(e "+ 1 - e "^|^ n+1 ,0), (104) 

where <£> is a real number, and that the energies of the supersymmetric partners H and H x 
are related by 

4 = e° n+1 . (105) 

Note that the operator x~ (respectively x + ) converts an eigenfunction of H (respectively 
Hi) into an eigenfunction of Hi (respectively H ) with the same energy. Thus, the opera- 
tors x~ and x + transfer the states from one spectrum to its partner spectrum. To introduce 
the ladder operators inside a given spectrum, we first consider the unitary operator U 
relating the states |^ n , 0) and |^ n , 1) through (cf ED- El) 

U := 1) (*n, 0| => |*„, 1) = U\* n , 0). (106) 

n 

Operators similar to U were already considered for continuous spectra II3T1 l33l and for 
discrete spectra 11321 l34l . Then, we define the operators (1321 - 11341 

a + ■= x + U a' := U ] x~. (107) 



19 



By using equations (11031) et (11041) . we obtain 

a-|*„,0) = v^e^-^-^I^^O) (108) 

a + |^„„0) = y^TTe-^-^^l^+^O). (109) 

Consequently, a + and a~ are creation and annihilation operators for the Hamiltonian H . 
Furthermore, it is easily seen that 

a+ a - = x + x~ = H . (110) 

Ladder operators for the Hamiltonian H\ can be introduced in a similar way. 

6.2 Physical realizations of the generalized oscillator algebra 

To simplify the notation, we set j^} := |\I/ n , 0). From equations (11081) et (11091 ), we get 

[a-,a + ]|tf n ) = (e n+1 -e n )|* n ). (Ill) 

The number operator N defined by 

N\* n )=n\V n ) (112) 

is in general (for an arbitrary quantum system) different from the product a + a~ . Let us 
consider the situation where the creation and annihilation operators satisfy the commuta- 
tion relation 

[a~,a + ] = aN + b, (113) 

a relation used in the study of the so-called polynomial Heisenberg algebra introduced in 
11361 . In other words, we assume that the energy gap & n +\ — e rt between two succussive 
levels is linear in n, i.e. 

e n +i - e n = an + b, (114) 

where a and b are two real parameters. We also assume that the eigenvalues of the operator 
aN + b are positive. With these choices, the algebra generated by the operators a + , a - 
and iV is identical to the generalized oscillator algebra A K modulo the replacements 

± a± la 
ar -» — = k — y — — (115) 

Vb 26 
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in equation (OQ). Thus, from equations (I108H110I) . we have 

H = a + a~ = ^aN(N - 1) + bN. (116) 
For a 7^ 0, the spectrum of H is non-linear and is given by 

H \^ n ) = e n \m n ) e n = -an(n-l)+bn. (117) 

Particular realizations of (11171) in terms of one-dimensional solvable potentials were pre- 
viously considered in [(9l[Tni[l2l[37l[38]|. Following the developments in lfT2l . we consider 
the following remarkable cases. 

(i) For (a = 0, b > 0), the spectrum of H is infinite-dimensional (n 6 N) and does 
not present degeneracies. 

(ii) For (a > 0, b > 0), the spectrum of H is infinite-dimensional (n & N) and does 
not present degeneracies. 

(iii) For (a < 0, b > 0), the spectrum of H is finite-dimensional with n — 0, 1, . . . , s— 
1 where 

7 Q 7 

s = — + -for -2- odd (118) 

a 2 a 

b b 

s = h 1 for - 2- even, (119) 

a a 

and all the states are nondegenerate. 

It is possible to find a realization of each of the three cases above in terms of exactly 
solvable dynamical systems in one dimension. We give below the corresponding potential 
v and transfer operators. 

(i) The case (a = 0, b = 1) corresponds to the harmonic oscillator (for which nGN) 
with 

v (x) = ^(x 2 - 1) (120) 



x ± = a ± = -4= ( =F-r- + x ) • ( 121 ) 



and 



a/2 V dx 

(For the harmonic oscillator, U reduces to the identity operator.) 

(ii) The case (a = 1, 2b = u + v + 1), with w > 1 and f > 1, corresponds to the 
Poschl-Teller potential (for which n G N) with 

— 1) f (f — l)l 1 



V {X) 



sin 2 | cos 2 | 
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u + vf (122) 



and 

i 1 [ d 1 / x x\~\ 
x = =F- — h — ( it cot — — f tan — ) . (123) 



d 1 / x x 

=F- — I — wcot t>tan — 

2V 2 2 



(iii) The case (a = — 1, 26 = 21 — 1), with / e N*, corresponds to the Morse potential 
(for which n = 0, 1, . . . , Z) with 

Uo(*) = ^ [e- 2 ^ - (2/ + l)e- a + Z 2 ] , (124) 



x ± = ^=(T^ + l-e- x ) . (125) 



and 



a/2 \ rfx 

6.3 Phase states and MUB for exactly solvable systems 

From equation (|53l , we can obtain the phase states for a general quantum system de- 
scribed by a truncated generalized oscillator algebra A K:S . We get 

1 5-1 

\m,<p) = -= Ve-^(g s )" m |vl> n ), (126) 

with s sufficiently large for the harmonic oscillator and the Poschl-Teller systems and 
s = I + 1 for the Morse system. Furthermore, equation (1871 ) provides with a mean to 
generate MUBs associated with the cases (i), (ii) and (iii) of subsection 6.2. 
On the other hand, the discrete phase state (l60l) reads here 

s-l 

<p) = C -= f e- fa »n&)"1*n>, ( 127 ) 

^ v^( n ) 

where the factor E{n) can be calculated in the different cases (i), (ii) and (iii). A simple 
calculation gives the following results in term of the Y function. 

(i) For the harmonic oscillator potential: 

E(n) = T(n+l). (128) 

(ii) For the Poschl-Teller potential: 

E{n) = r '" + W" + " + | " +1 » . (129) 
2 n T(u + v + 1) 



(iii) For the Morse potential: 



= n.;TO, (130) 

v ; 2 n T(2l-n) 
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It should be mentioned that the discrete phase states given by (|127l) differ from the 
coherent states for exactly sovable potentials derived in ifTOl l34l [36l l38l l39l from super- 
symmetric quantum mechanics techniques. The noticeable difference comes from the fact 
that the states (11271) are temporally stable and are labeled by an integer instead of a contin- 
uous complex variable as in the coherent states derived in [[TOl [34l [36l [381 [39J]. The states 
(11271 ) are eigenstates of the operator (1551 ) whereas the coherent states in [fT0l[34l[36l[38l[39ll 
are obtained from the three standard definitions (involving annihilation operator, displace- 
ment operator, and uncertainty relation). 

7 Concluding remarks 

The starting point of this article is based on the definition of a generalized oscillator 
algebra A K . This algebra is interesting in two respects. First, it describes in an unified 
way some exactly solvable one-dimensional systems having a nonlinear spectrum (for 
k 7^ 0) or a linear spectrum (for k = 0). As typical examples, these quantum systems 
correspond to the Poschl-Teller potential (for k > 0), the Morse potential (for k < 0) 
and the infinite square well potential (for k = 1/3) in addition to the harmonic oscillator 
potential (for k = 0). Second, the algebra A K can take into account some nonlinear effects 
that may occur in the quantum description of quantized modes of the electromagnetic field 

(cf. ed. 

In connection with the algebra A K , the present work adresses three problems: the 
construction of a phase operator, the determination of its temporally stable eigenstates (the 
so-called phase states) and the derivation of MUBs from the obtained phase states. This 
is the first time that a connection between MUBs and dynamical systems is established. 
In this regard, the character "temporally stable" of the eigenstates of the phase operator is 
essential for the derivation of MUBs. The main results of this paper are as follows. 

For the case k > (which corresponds to an infinite representation of A K ), the phase 
operator is not unitary. We note in passing that the corresponding phase states are similar 
to those derived in ifTTl except that our states are temporally stable. However for k > 0, 
by making a (a la Pegg and Barnett) truncation, which gives rise to a truncated generalized 
oscillator algebra A KjS , we can define a unitary phase operator whose eigenstates lead to 
MUBs. 

For the case n < (which corresponds to a finite representation of A K ), it is possible to 
construct a unitary phase operator whose eigenstates are temporally stable. MUBs can be 
derived as a subset of these states. For k < 0, the consideration of a truncated generalized 
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oscillator algebra A KyS is nevertheless necessary in order to establish a connection with 
the Morse system and to derive associated MUBs. 

As a conclusion, in both cases (k > and k < 0), the truncation procedure makes it 
possible to define a unitary phase operator for exactly solvable systems and to generate 
temporally stable phase states from which MUBs can be derived. 

Another result of this paper concerns a new type of phase states. These temporally 
stable phase states, namely the states (|60l) . are associated with the truncated algebra A K , S - 
They are eigenstates of an operator defined in the enveloping algebra of A K jS and con- 
stitute discrete analogs of the coherent states derived in [8]. More generally, this result 
shows that it is possible, for a finite spectrum, to derive new phase states similar to the 
coherent states of (8]| constructed, for an infinite spectrum, as eigenstates of an annihila- 
tion operator. The key of the derivation of the new states (for a finite spectrum) is to add 
a power of the creation operator to the annihilation operator. 

To close this paper, let us mention that the concept of MUBs was recently extended 
to infinite-dimensional Hilbert spaces iflTTl . In this vein, it is hoped that the temporally 
stable phase states derived in this work for the infinite-dimensional case could serve as a 
hint for deriving MUBs for continuous variables, a difficult challenge. 
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